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I. INTRODUCTION 

In the study of gravitational collapse, we always encounter with two types of singularities 
— shell focusing singularity and shell crossing singularity. In Tolman-Bondi-Lemaitrc 
(TBL) dust model, these two kinds of singularities will corresponds to R = and R' = 
respectively. A shell focusing singularity (i.e., R — 0) on a shell of dust occurs when it 
collapses at or starts expanding from the centre of matter distribution. The instant at 
which a shell at the radial co-ordinate r will reach the centre of the matter distribution 
should be a function of r. So different shells of dust arrive at the centre at different times 
and there is always a possibility that any two shells of dust cross each other at a finite 
radius in course of their collapse. In this situation the comoving system breaks down, 
both the matter density and kretchmann scalar diverge [1,2] and one encounters the shell 
crossing singularity. If one treats it as the boundary surface then the region beyond it is 
unacceptable since it has negative density. It is therefore, of interest to find conditions 
which guarantee that no shell crossing will occur. 

In TBL model, shell crossing singularity has been studied by several authors [1-9]. Also 
Goncalves [7] studied the occurrence of shell crossing in spherical weakly charged dust col- 
lapse in the presence of a non- vanishing cosmological constant. The positive cosmological 
constant model conceively prevent the occurrence of shell crossing thereby allowing at least 
in principle for a singularity free 'bounce' model. Nolan [8] derive global weak solutions 
of Einstein's equations for spherically symmetric dust-filled space-times which admit shell 
crossing singularities. Recently, Hellaby et al [9] investigate the anisotropic generalization 
of the wormhole topology in the Szekeres model. In this work, we have studied the shell 
crossing singularity in Szekeres model of the space-time both from physical and geometri- 
cal point of view. In section II, we derive the basic equations in shell focusing and shell 
crossing singularities. We study the physical conditions and geometrical features of shell 
crossing singularities in sections III and IV respectively. Finally the paper ends with a short 
discussion in section V. 
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II. BASIC EQUATIONS IN SHELL FOCUSING AND SHELL CROSSING 

SINGULARITIES 

Recently, we have presented dust solutions for (n + 2)-dimensional Szekeres' space-time 
model with metric ansatz [10] 

n 

ds 2 = dt 2 -e 2a dr 2 -e 2 ^^ W 

i=l 

where a and (3 are functions of all the (n + 2) space-time co-ordinates. If we assume that 
f3' (= || ) 7^ 0, then the explicit form of the metric coefficients are 



and 



where 



e? = R{t,r) e^i.-^n) (2 ) 



e" - R ' + RU ' (3) 



e-^iW^^+^B^+CM (4) 



»=i »=i 
and i? satisfied the differential equation 

Here A is the cosmological constant, /(r) and F(r) are arbitrary functions of r alone; and 
the other arbitrary functions, namely A{r), Bi(r)'s and C(r) in equation (4) are algebraically 
related by the relation 



Y J B 2 -AAC=-l. (6) 

i=l 

It is to be noted that the r-dependence of these arbitrary functions A, Bi and C play an 
important role in characterizing the geometry of the (n + l)-dimensional space. In fact, the 
choice A(r) = C(r) — | and Bi{r) = (V i = 1, 2, n) reduce the space-time metric (1) to 
the usual spherically symmetric TBL form 



ds 2 = dt 2 - ±- dr 2 - R 2 dn 2 n . 
1 + f(r) 

by the co-ordinate transformation 

X\ = Sin9 n Sin9 n _i Sin9 2 Cot^9i 

x 2 = Cos0 n Sin9 n _i Sin9 2 Cotj9i 

xs = Cos6 n -iSin6 n -2 Sin9 2 Cot\9\ 



x n-i = Cos9^Sin9 2 Cot^9i 
x n = Cos9 2 Cot^9i 
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Hence in the subsequent discussion we shall restrict ourselves to the quasi-spherical space- 
time which is characterized by the r dependence of the function v (i.e., v' ^ 0). The 
expression for energy density due to dust matter using the Einstein equations is 

N n F' + (n + l)iV 
P (t,r, Xl ,...,x n ) = - Rn{RI + Rvl) ■ (8) 

The space-time singularity will occur when either (i) R = i.e., {3 — — oo or (ii) 
R' + Rv' = i.e., a = — oo. The standard terminology for spherical collapse suggests that 
the first case corresponds to shell focusing singularity while in the second case we have a 
shell-crossing singularity. In the following we shall discuss the situations for shell crossing 
singularity. 

Suppose the collapse develops at the initial hypersurface t = U where we assume R(ti, r) 
to be a monotonically increasing function of r. So, without any loss of generality, we can 
label the dust shells by the choice R(t i7 r) = r. Hence the expression for the initial density 
distribution is given by 

Pi{r,xi, ...,x n ) = p{U,r,xi, ...,x n ) = - — — r — (9) 

2 r"(l + rv') 

If we consider that the collapsing process starts from a regular initial hypersurface then the 
function pi must be non-singular (and also positive from physical point of view). Moreover 
the local flatness property of the space-time near r — demands f(r) — > as r — > 0. Then 
in order to R 2 to be bounded as r — > we must have F(r) ~ 0(r m ) where m > n — 1 (see 
eq. (5)). On the other hand, for small r, Pi(r) ~ f F +( '"t 1 ' > ' F ' y and consequently, for regular 
pi(r) near r = 0, we must have F(r) ~ 0(r n+1 ) and v' ~ O(^). Hence, starting from a 
regular initial hypersurface, we can express F(r) and Pi(r) as a power series near r — as 
[11] 



F{r)=Y J Fjr n+i+l (10) 

and 

oo 

Pi(r)=J2pi ri - ( n ) 

As i/ appears in the expression for the density as well as in the metric coefficient, so we 
can write [11] 



v'{r) = J2 "J (I 2 ) 
j=-i 

where v_, > —1. 



Now, using these series expansions in equation (9) we have the following relations between 
the coefficients, 

n(n + 1) _ n ( , 1 \ _ 



and so on. 



Now in order to form a singularity from the gravitational collapse of dust, all portions of 
the dust cloud should collapse i.e., R < 0. Let us denote by t s f(r) and t sc (r) as the time 
for shell-focusing and shell-crossing singularities occurring at radial coordinate r. Hence we 
have the relations 

R(t sf ,r) = (13) 

and 

R'(t sc ,r) + R(t sc ,r)u'(r,x 1 ,X2, ...,x n ) = 0. (14) 
Note that 'i sc ' may also depend on xi, X2, x n . 



III. PHYSICAL CONDITIONS FOR SHELL CROSSING SINGULARITY 

We shall now make a comparative study of shell focusing and shell crossing singularity 
time and find conditions in favour (or against) of formation of shell crossing singularity for 
the following different choices: 

(i) /(r) = 0, A = 0: 

In this case equation (5) can be integrated to give 

R ^l^ rH l_n + l^ ¥F){t _ u) (15) 

So R(t s f,r) = results 



tsf(r) =U + ———=== r~ (16) 
(n + l)^/F(r) 

Now to avoid the shell crossing singularity either all shells will collapse at the same time 
(i.e., t s f is independent of r) or larger shell will collapse at late time (i.e., t s f(r) is a monotone 
increasing function of r). These two conditions can be combined as 

t' sf (r)>0 

or equivalently from equation (16) 



F(r) 

Now combining equations (14) and (15) we have 



F'(r) < n + 1 
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t sc (r) - t 3f (r) = 



(n + l)VFM{^g + ( n +iy) 



(18) 



But if it is so happen that R' + Rv' = is a regular extremum for (3, then we must 
have finite p. This implies from equation (8) that F' + (n + l)Fu' = 0. Hence, if there 
is no shell crossing singularity corresponding to equation (14) we must have two possibilities: 



F'(r) n 4- 1 F'(r) 

either (a) tL^l + (n+ l)z/ = and -— - -±1 = Q 



F(r) 



F(r) 



(19) 



or 



(b) m + {n + iy = and IL±1-S!2>„ 



F(r) 



F(r) 



(20) 



For the first choice i s / is constant, so all shells collapse simultaneously while for the 
second choice t s f is a monotonic increasing function of r and there is an infinite time 
difference between the occurrence of both type of singularities. 

The value of R at t = t sc (r) is 



n±i f F'( r) _ n+1 -j 
+ („+!),/} 



(21) 



Therefore as an complementary event, the conditions for occurrence of shell crossing 
singularity are R' + Rv' = 0, p = oo, R < 0, i? > 0, i' s/ < 0. 



As t' s j < implies 



so R > demands 



F'(r) n + 1 
> r 



(22) 



Hence we have 



F(r) 



+ (n+ l)v' > 



F(r) - r l and e 1 ' - r p 
for shell crossing singularity with I > (n + 1) and (n + l)f> > — I. 



(23) 



(24) 



(ii) /(r) = 0,A^0: 

This choice will give the solution to equation (5) as [11] 



t = U + 



2n 



(n + l)A 



Smk " (\/ n(n 2 +l)F(r) ) " (^ n(n^(r) 



(25) 



G 



At the shell focusing time t s f(r), R = 0, hence we have 



tsf(r) 



2n 



(n + l)A 



Sink 1 



/ 2Ar n+1 
n{n+l)F(r) 



(26) 



Since to avoid the shell crossing singularity, t s f(r) should be an increasing (or a constant) 
function of time i.e., t' s ^(r) > 0, which implies the same condition (17) as in case (i) and is 
independent of A (whether zero or not). Further, substitution of equation (25) in equation 
(14) will give 



t sc {r) - t s f(r) = tank 1 



.{f^ + (» + lK} 



(27) 



But if there is no shell crossing singularity corresponding to R! + Rv' = (then it will 
correspond to an extremum of (3) then p must be finite. This will be possible only when 
F' + (n + l)Fv' = 0. But from cq.(27) it is permissible only when t' s j = i.e., t s j is 
independent of r. 



(hi) f(r) 7^ 0, A = 0, R(t i} r) — (time symmetry): 

In this case explicit solution is possible only for five dimension (i.e., for n = 3) and the 
result as 



R 



2 _ 2 F(r) 



(t - Uf 



(28) 



But the shell focusing condition R(t s f,r) — gives 



t s f — U + 



So t' s f > will give 



F 4 
— < - 

F ~ r 



Here the time difference between the two types of singularities is 

1 



tsc tsf 



2F x " r 



(29) 



(30) 



The r.h.s. of equation (30) always positive by the inequality (29). 



IV. GEOMETRICAL FEATURES OF SHELL CROSSING SINGULARITY 

Now wc shall discuss the shell crossing singularity from geometrical point of view. We 
note that R' + Rv' (related to e a ) is always > 0. The equality sign corresponds to the shell 
crossing singularity. In fact geometrically, a shell crossing singularity (if it exists) is the 
locus of zeros of the function R' + Rv' (i.e., a = — oo). Now writing explicitly the function 
R' + Rv' using the solution (4) for v we have 



R' + Rv' = e" 



(R'A - RA') J2 ^ + J2( R ' Bi ~ RB 'i) Xi + ( CR> ~ C ' R ^ 



i=i 



i=i 



(31) 
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We see that (a detailed analysis is given in the Appendix) there will be no shell crossing 
singularity i.e., R' + Rv' will be positive definite if 



n/2 " 

-^>5>?-4A'C' = V>(r)(say) (32) 

i— 1 

Note that R' + Rv' will also be positive for -^j- = V'( r )i provided Xi ^ ^w'a-ra 1 ) = 
xoi , i = 1,2, ...,n. Thus shell crossing is a single point (xoi, £02, ...,xq u ) in the constant 
(t, r)-hypersurface (n dimensional). In other worlds, it is a curve in the ^-constant (n+ 1)-D 
hypersurface and a 2 surface in (n + 2)-D space-time. 

When < ip(r) then shell crossing singularity lies on n-hypersphere in the n-dimensional 
Xi's plane. This hypersphere has centre (xqi, xq2, ...,Xo n ) and radius 



2(i?M - JL4') 1 j 

In the above we have assumed a = — RA') to be positive. However, if a < and 

-^2" < V'( r ') then also shell crossing singularity is possible and it lies on an n-hypersphcre 
having same centre (xoi, X02, Xo n ) an d radius r c . Further, there will be no shell crossing 
singularity if the variables Xj's lie inside the above n-hypersphere. 

The above hypersphere is different from the hypersphere with v' = i.e., 



A'(r)J2 x i+J2 B 'i^ Xi + C '^ =0 ( 34 ) 

i=l i=l 

So the shell crossing set intersects with the surface of constant r and t along the line 
(curve) ^ = — i/=constant. 

Now for positive density we note that F' + (n + l)Fv' and R 1 + Rv 1 must have the same 
sign. We now consider the case where both are positive (when both are negative, we just 
reverse the inequalities). When both are zero then it can happen for a particular value of 
Xi's (i = 1,2, ...,n) if r^p\Tp = 7^" = — v> i which can not hold for all time. This is possible 
for all x t if F' = R' = v' = 0. This implies that at some r, F' = f = A' = C = B[ = 
(i = 1,2, ...,n). Hence we choose 

F' R' 

> — v' and — > -v 1 (35) 



(n+l)F R 
Also from the solution (4) we have 

_ , _ A'(r) Y.U A + E?=i B[(r)x t + C'(r) 
A^Y.U^ + Y.UB^x. + Cir) 

Now writing in a quadratic equation in x\ we have for real Xi, 



v n + 2(A> + Ay') x k + + B k y') < ]T B? AA'C 

{ fc=2 k=1 J i=l 

So 

n 

v a \ max = Y B ?- AA ' C ' 



i=l 
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Hence from (35) we have 



F' 



> 



{ n+l)F-\^ 



B? - iA'C , V 



(36) 



which implies F' > 0, V r. Now for R' + Rv' > 0, we shall study the three possible 
choices separately. 

(i) /(r) = 0,A = 0: 

Here the solution for R can be written as 

R^ = i ^VW){t-a{r)) 

Soast^a, R^R' + R^v' -> -y/F(r) a'(r) and as t -» oo, ^ + 1/' -» ^jjp + v' ■ 
Hence for i?' + Rv' > we must have a' < and (w ^ )F > VE^i B f = 4 ^'C" • 



(ii) /(r)=0,A = 0: 

The solution for R can be written as 



R^ = 5in * 
2A 



(n + 1)A 
2^ 



(i-a(r) 



In this case as t — > a, R™? 1 R' + R™^ v' — > —y/F(r) a'(r) and as i ^ oo, 



(n+l)F 



o'(r) + i/'. Thus for i?' + > we must have 

a > {r) < o and ^ J^Ml a , {r) > VEr=i^-4A'C . 



(iii) /(r) 7^ 0, A = 0, i?(t*,r) =0, n = 3: 
Here the solution for R is 



# = r 2 - F(r) - < 



The limiting value of ^ + z/ as i — > oo will be 2 , 
should have |£ + ± > VE"=i B f - ^ A ' C ' ■ 



F' 



: . Hence for R' + Rv' > we 



DISCUSSION 



In the last two sections a details study of shell crossing singularity has been done for dust 
model with or without cosmological constant for Szekeres model of (n + 2)-D space-time. 
The physical conditions for shell crossing singularity are presented in section III. These 
conditions however do not depend on A (whether zero or not) and the form of the conditions 
are identical for the three cases presented there. For geometrical conditions the locus of 
shell crossing depends on the discriminant of the co-ordinate variables x^s (i — 1,2, ...,n). 
If both ^ and jj^ijp are greater than yE B'^ ~~ 4A'C then there will be no shell crossing 
singularity even if R' + Rv' = 0. Here p is finite and R' + Rv' — will correspond to 
a real extrema for /3. On the other hand if ^ = \/E B'^ ~ 4A'C then shell crossing 
singularity is a 2-surfacc in (n + 2)-dimensional space-time. For ^ < B'^ — AA'C, 
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the shell crossing set lie on a n-hypersphere and it intersects with constant (t, r) along the 
curve ^ = — v' = constant. For future work it will be interesting to study in details the 
possibility of shell crossing singularity with pressure. 



Appendix: A detailed study of a quadratic expression: 

Consider a general quadratic expression in n variables 



n n 



z = a Xi + b i X i + c ( 37 ) 
i=l i=l 

For the present problem (given in equation (31)) we have 

a = R'A - RA' , bi = R'Bi — RB[ , c= CR' - C"R (38) 
Equation (37) can be rewritten as 



n / , \ 2 , n 

i=l v 7 i=l 



(39) 



By substitution from (38) we obtain 

d = MR 2 + NR' 2 + LRR' (40) 

with 

n n n 

M = AA'C' -J2 B 'i 2 ' N = AAC ~ H B i ' L = -4(A'C + AC') + 2Y B i B 'i 

i=l i=l i=l 

Using equation (6) one sees that N = 1 and L = and d simplifies to 

d = MR 2 + R' 2 
We shall now discuss the sign of z for the following cases: 

(i) a > 0, d > : It is clear from equation (39) that z is positive definite for all values 

of Xi's. 

(ii) a > 0, d — : In this case z > 0. The equality sign occurs for a particular value of 
the variables Xi's namely xi = — |^ , i = 1, 2, n. 

(m) a > 0, < : Here z has indefinite sign. In particular, z will zero when the 
variables will lie on a hypersphere having centre (— — |j) an d radius 

If a < we have indefinite sign of z for d < and for d > 0, z will be negative definite 
(which is not possible in the present paper). We note that if z = (with d < 0) then as 

above the variables Xi's lie on the hypersphere having the same centre and radius is -^-5^. 
But if the variables lie inside the above hypersphere then z will have positive value. 
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